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Abstract 

So called pair copula constructions (PCCs), specifying multivariate distributions 
only in terms of bivariate building blocks (pair copulas), constitute a flexible class of 
dependence models. To keep them tractable for inference and model selection, the 
simplifying assumption that copulas of conditional distributions do not depend on the 
values of the variables which they are conditioned on is popular. In this paper, we 
show for which classes of distributions such a simplification is applicable, significantly 



extending the discussion of |Hobaek Haff et al. (2010) In particular, we show that the 



only Archimedean copula in dimension d > 4 which is of the simplified type is that 
based on the gamma Laplace transform or its extension, while the Student-t copula is 
the only one arising from a scale mixture of Normals. Further, we illustrate how PCCs 
can be adapted for situations where conditional copulas depend on values which are 
conditioned on. 

Keywords: archimedean copula; elliptical copula; pair copula construction; condi- 
tional distribution 



1 Introduction 



Growing capabilities to store large data sets and increasing computing power for 
their computational analysis substantially increased the demand to develop statisti- 
cal methodology for many dependent variables over the last years. The properties of 
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2 Pair copula constructions 



high-dimensional distributions in general, however, remain hard to understand intu- 
itively. Given that, the early focus was on the multivariate Normal distribution, which 
allows to describe the whole dependence structure by specifying bivariate correlations. 
These linear dependencies between just two variables remain easy to understand and 



communicate. But more recent developments, in particular in econometrics (c.f. Lon- 



gin and Solnik (2001) , Ang and Chen (2002) ), showed that many data sets in higher 



dimensions exhibit more complex dependence structures. 

To deal with these challenging structures and keeping the model understandable at 
the same time, a popular technique is to sequentially decompose the joint distribu- 
tions into bivariate building blocks by conditioning. These are then modeled using 
bivariate copulas of which many parametric families are well studied, see for exam- 



ple the monographs by Joe (1997) and Nelsen (2006) Distributions arising from 



such a decomposition are called pair copula constructions (PCCs). The question for 
which classes of multivariate models and under which assumptions such a sequential 
conditioning is possible, however, remained unsolved. The only publication in this 
direction is Hobaek Haff et al. (2010) , providing first illustrative examples. 

In this paper, we fill in this gap by classifying multivariate distributions in terms of 
the dependence properties of their conditional distributions (see Figure [TJ. Special 
attention will be paid to the simplifying assumption that the copulas corresponding 
to conditional distributions are constant irrespective of the values of variables that 
they are conditioned on. This assumption is usually made to keep model selection 
and inference tractable, and we will show for which multivariate distributions it holds. 
Furthermore, we discuss how constructions based on simplified PCCs can be extended 
when the simplifying assumption is not applicable. 

The remainder is structured as follows. Section[2]introduces PCCs. In Section [3] and |4| 
we demonstrate which Archimedean and elliptical copulas are simplified PCCs. Sec- 
tion [5] provides the remaining examples for the classification in Figure [T] and considers 
the increased flexibility gained by weakening the simplifying assumption. Section [6] 
concludes with an outlook to areas of future research. 



2 Pair copula constructions 

A copula is defined as the cumulative distribution function (cdf) of a multivariate 
distribution with uniform(0, 1) margins. These have become of great importance 



in multidimensional dependence modeling, since the theorem of Sklar (1959) allows 



to separate the specification of a d-dimensional distribution F into specifying the 
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Figure 1 Overview of the PCC classes which are discussed including their relations 
and examples. 

univariate marginal distribution functions F^, i = l,...,d, and a copula C, i.e., to 
find C such that 

Fi !d (a;i, ...,x d )= C 1 . A {F 1 (x x ), . . . , F d (x d )). 



Here, we use the abbreviation i : j := + for i < j. In particular 

bivariate copula families are well studied, and every multivariate distributions can be 
constructed from marginal distributions and two-dimensional copulas by sequential 
conditioning, a so-called pair-copula construction. While different principles for PCCs 



such as regular vines (R-vines) (Bedford and Cooke (2001, 2002)) or non-Gaussian 



directed acyclic graphs (DAGs) (Bauer et al. 2011) have been investigated over the 



years, the underlying basic idea is always that of Joe (1996) 



Let us consider three random variables (rvs) Xi, X 2 and X 3 with corresponding 
cumulative distribution functions (cdfs) Fi,F 2 and F 3 , then we can first construct a 
joint distribution of (Xi,X 2 ) and of (X3, X2) by assigning copula functions C\ 2 and 
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2 Pair copula constructions 



C32, respectively: 

F 12 (xi,x 2 ) :=C 12 (F 1 (x 1 ),F 2 (x 2 )) 
F 32 (x 3 ,x 2 ) := C 32 {Fx (x 3 ), F 2 {x 2 )) 

These bivariate distributions can be combined to a three-dimensional distribution by 
assigning a conditional copula C1312, 



F 123 {xi, x 2 , x 3 ) := / Ci 3 \ 2 (Fi\ 2 (xx\x 2 ),F 3 \ 2 (x 3 \x 2 );x 2 ) dF 2 (x 2 ). (2) 
J —00 

For every cdf -F123, Sklar's theorem implies the existence of a bivariate copula that 
couples Fi\ 2 (-\x 2 ) and -^312 (' I ^2) for every x 2 . The generalization of this principle to 
arbitrary dimensions is obvious and the conditional copulas corresponding to such a 



construction can be determined for every multivariate distribution (c.f. Patton (2006) 
for a discussion of conditional copulas). 

To make PCCs tractable for inference and model selection however, further assump- 
tions regarding the conditional copulas have to be made. 

We will call a multivariate distribution an absolutely continuous and parametric 

PCC if all bivariate copula families occurring in the construction have densities with 
a parameter vector. In this case, it will be possible to express the likelihood function 
of the joint copula model in terms of the bivariate densities and to make inference 
about the finite-dimensional parameter based on this expression. One example for 
this would be the decomposition of a <i-dimensional density / with marginal densities 
fk, 1 < k < d, as a C-vine copula (see Aas et al. (2009)] ): with univariate parame- 



ters 6k and copula parameters and 6 = (6*i, . . . , 6d', — 1, ■ ■ ■ , d — 1, i 

l,...,d-j), 



f 1:d {xi, , x d ; 0) = f k {x k ; 6k) ■ 

k=l 

d-l d-j 

II II C JJ+i|l=0-l) ( jF i|l:(i-l)( a? J ^+^11:0-1) (^i+i I X l : (j- 1 ) ) 5 ^'j+O ■ 



j=l i=l 

The above multivariate density is valid for a simplified PCC because all bivari- 
ate copulas occurring in the constructions do not depend on the variables that are 
conditioned on, e.g., 

C\ 3 \ 2 {-, x 2 ) = Ci 3 | 2 (-, •)• 

This simplifying assumption reduces the specification of a PCC to choosing bivari- 
ate copula families and their parameters. This means that the potentially complex 
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3 Archimedean copulas 



dependence between variables that are conditioned on and the copula functions can 
be neglected, and model selection can be performed more easily. Also other tech- 
niques for PCCs such as for R-vines (stepwise estimation in Hobaek Haff (2012) and 
Difimann et al. (2011) , inference as in Aas et al. (2009) or Bayesian techniques as in 



Min and Czado (2010) ) extend to the non-simplified case only if the dependence on 



the variables that are conditioned on is of a simple form (we discuss this in Section 
[5l see in particular Example 5.1). 



3 Archimedean copulas 



In this section, we characterize the Archimedean copulas that are simplified PCCs. 
A G?-dimensional Archimedean copula is given by 

d 



C 1:d (ui, ...,u d ) 



i=i 

for an Archimedean generator function ip G C^- Here, is the class of functions 
if : [0,oo) 1 — y [0, 1], which are strictly decreasing on [0, inf {s, (p(s) = 0}] (with ip(0) = 
1, (p(oo) = 0) and differentiable on [0, 00) up to order d — 2, such that (— > 
0, j = 1, ...,d — 2, and where further (— l) d ~ 2 Lp^ d ~ 2 "> is non-increasing and convex 



(3) 



(McNeil and Neslehova (2009) , Joe (1997) , Miiller and Scarsini (2005) ). These copulas 



have the appealing property that the conditional distribution function is given in a 
simple analytical form which facilitates the study of their properties. Given a Gamma 
random variable with shape parameter p, rate parameter b, density 



fr(x;b,p) 



r(p) 



mean p/b and variance p/b 2 , the Laplace transform (LT) is 



°f r {x;l/0,l/9)dx 



-i/e 



6>0. 



The Archimedean copula corresponding to this LT as generator function is called 
MTCjQ copula. From the form of the conditional distributions derived in Takahasi 



1 This copula is also called Clayton copula due to its appearance in ( Clayton 1978 ) . It is the copula 



oTThT^rmrH^Mate*F^areTo"T^^ 



(|Takahasi 1965[). It was first mentioned as a multivariate copula in|Cook and Johnson (1981)|and 
as a bivariate copula in|Kimeldorf and Sampson (1975) The extension to negative dependence 
was given in Genest and MacKay (1986)|for the bivariate case and in Joe (1997, pp. 157-158) 



for the multivariate case. Since many properties were discovered studying the corresponding 



distribution function and Cook and Johnson (1981) mentioned its general form we call it MTCJ 
copula — 
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(1965) it is obvious that the copulas corresponding to bivariate conditional margins 
of the MTCJ copula are again MTCJ copulas. The above gamma LT family extends 
to Archimedean generators (for the generalized MTCJ family) 

M^0) = (l + O-s)- 1/e , 0>-^lp (:r) + = max{0,a;}. (4) 



The conditional generator (c.f. Mesfioui and Quessy (2008) ) when conditioning on m 
variables becomes 

g Q N -(me+l)/0 



1 + ^nJ + • (5) 

For 9 > l/(d — 1), 6/(m6 + 1) > l/(d— 1 — m), so the parameter range is consistent, 
and the (generalized) MTCJ copula is a simplified PCC where all bivariate building 
blocks are MTCJ with corresponding choice of parameters. In fact, under a weak 
regularity assumption, it is the only multivariate Archimedean copula that constitutes 
a simplified PCC. 

Theorem 3.1 

A d- dimensional Archimedean copula with a generator that is twice continuously dif- 
ferentiable on the set where it is positive is a simplified PCC if and only if its generator 
is in the family Q). 



Proof 

Appendix 



A 



□ 



Note that the differentiability condition for the generator is satisfied for all Archimedean 
copulas if d > 4. This adds a further aspect making the MTCJ copula unique^] among 
Archimedean copulas. Using different parameters than obtained through Equation 
(J5|, PCCs with MTCJ copulas with freely chosen parameters as building blocks nat- 
urally generalize the MTCJ copula to what we call an extended MTCJ copula. 



4 Elliptical copulas 

In this section, we characterize the elliptical copulas that have all conditional distribu- 
tions in location-scale families. We show that not all elliptical copulas are simplified 
PCCs and characterize the scale mixtures of Normals which are of the simplified type. 
By referring to a d- dimensional elliptical copula we mean a copula arising from an 



2 The MTCJ copula also is the only Archimedean copula invariant under truncation, in the s ense 
that for a rv U - C, C also is the copula of U|U < a, a e [0, l] dim(u) (jAhmadi Javid 20091. 
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elliptical distribution such as the multivariate Normal or Student-t distribution. A 
multivariate distribution is elliptical if its characteristic function has the form 

9 x (t; A i,I!) = *(t'£t)e*''*, 

for \& : Rq~ 4l, G M. d , and positive definite S G R dxd . If the distribution has a 
density, this implies that it is given by 



/x(x) = |S|- 1 /2^ (x _ / ,y s - 

for a generator function g : ^ ^o~; which can be uniquely determined from see 



Cambanis et al. (1981) For the two examples we mentioned, the generator functions 



have the form 

and lead to simplified PCCs. 
Theorem 4-1 

The multivariate Gaussian distribution and the multivariate Student-t distribution are 
PCCs of the simplified form. 

For the Gaussian distribution this is quite obvious since all conditional distributions 
are again Gaussian and do depend on the values of the variables that they are con- 
ditioned on only through their mean, i.e., they all are in the same location family. 
For the Student-t distribution the covariances also depend on these values, but only 
through a scaling factor such that all conditional distributions remain in the same 
location-scale family. Since changes of location and scale do not affect the copula of 
a multivariate distribution, this implies that the Student-t distribution is a simplified 
PCC. For clarification, we derive the explicit form of the copulas corresponding to 



bivariate conditional distributions in Appendix B.l Just as for the MTCJ copula, 
changing the degrees of freedom in the bivariate building blocks of the PCC leads 
to a natural extension of the multivariate Student-t distribution in which different 
bivariate conditional margins can have different degrees of freedom. 

While the copulas corresponding to these distributions are the most popular examples 
they also have a unique position within the class of elliptical distributions as the 
following theorems show. 
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4 Elliptical copulas 



Theorem J±.2 

Let us assume that the generator function g(-) of the density of a d- dimensional ellip- 
tical distribution is differentiable. Then, the conditional distributions remain within 
the same location- scale family for all values of the variables that are conditioned on, 
if and only if 

(a) the support of g isM. and the distribution is the multivariate Student-t (or Pear- 
son type VII) distribution or in its limiting case the multivariate Normal distri- 
bution or 

(b) g has compact support and g{t;Q = (1 — t)+, for £ > 1 up to rescaling (the 
Pearson type II distribution). 



Proof 

Appendix 



B.2. □ 



Here, case (b) can be seen as an analogue to the extension of the MTCJ to negative 
dependence. 

The proof that the Student-t distribution is a simplified PCC relied on the fact that in 
this case, conditioning only affects the location and scale of the distribution. Theorem 



4.2 now shows that the t-distribution is the only elliptical distribution where this proof 
strategy is successful. Just as the MTCJ can be constructed using a Gamma mixture, 
also the Student-t distribution arises from the multivariate Normal distribution using 
a Gamma mixture for the square of the inverse scale parameter. For the distribution 
function of the MTCJ copula, we have 

poo d 

C M TCj{^i:d\ 0) = JT [G MT cj(uh 6)V /r(«; 1, l/9)da, 
Jo i=i 

where Gmtgj{'\ @) — ex P {~ {fMTCj}^ 1 (s &)) is a cdf on [0, 1], c.f. Joe (1997, p. 86) 
while 

POO 

/ M (x 1:d ; R, v) = (2n)- d ^\R\-^ 2 / w d/2 exp{-iu«'ir 1 x}/ r («;; v/2, 2/v) dw, 

Jo 

holds for the density of a multivariate Student-t distribution. When only considering 
scale mixtures of Normals, we obtain the following theorem. 

Theorem 4.3 

Consider a d- dimensional scale mixture of Normals which is a simplified PCC 
a) in d > 4 or 
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b) for all correlation matrices S 7 



then it is a simplified PCC if and only if the mixing distribution is the Gamma dis- 
tribution. 



Proof 

Appendix 



B.3 



□ 



This contradicts the claim made in Example 4.1 of Hobaek Haff (2012) that all ellip- 
tical distributions are simplified PCCs. 

Although this shows that the intersection between simplified PCCs and elliptical dis- 
tributions is limited, it turns out that the deviation from the simplifying assumption 
can be ignored in many applications. In fact, the values of Kendall's r corresponding 
to bivariate conditional distributions are independent of the variables that are condi- 



tioned. Cambanis et al. (1981) showed that the conditional correlation coefficient is 



equal to the partial correlation for elliptical distributions, and Lindskog et al. (2003) 



demonstrated that the relationship between the correlation coefficient p and Kendall's 
t for the Normal distribution 



7T 



-arcsin(p), 



holds for all atom-free elliptical distributions. 



5 Non continuous and non-simplified PCCs 

In order to complete the list of examples for the classification in Figure [TJ we need 
to find distributions which are not absolutely continuous but simplified PCCs and 
distributions which are neither absolutely continuous nor simplified PCCs. For the 
first case, let us consider a PCC with bivariate Cuadras-Auge (which are a special case 
of the bivariate Marshall-Olkin (MO) copula, where the distribution is exchangeable) 
copulas as building blocks. Exempli gratia, 

612(^1,112) = min(Mi,M 2 ) max(Mi,M 2 ) 1 ~ ai2 , CK12 € [0,1], 



c.f. Mai and Scherer (2012, Chapter 1) For the joint distribution of (Ui, £/ 2 , U3), -F123, 



this implies that 



/ mm(F ai2 (ui\v 2 ),F a:i2 (u 3 \v2)) max(F ai2 (ui\v 2 ) , F a32 (u 3 [i^)) 1 Ql3|2 cfo 2 
J 
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where the conditional cdf is given by (Mai and Scherer 2012, Example 1.8) 

(u(l - a)v~ a u < v, 
F a (u\v) = < 

\u a U > V. 

Note that here F123 is a copula. A scatterplot of the bivariate 13 margin of this 
distribution is shown in the left panel of Figure [2] for illustration. 

While this copula is of simplified form by construction, the multivariate MO copula 



c.f. (Mai and Scherer 2012, Chapter 3)) 



C MO {ui:d)= II ™m{uf J:k * jXj \, Aj>0, £ Aj>0 Vfc = l,...,d, 

in contrast is neither a simplified nor an absolutely continuous PCC. We derive the 
copula of bivariate conditional distributions for the trivariate case in Appendix O 
the cdf of the conditional copula of Ui, U2\Us = u% is plotted in the right hand panel 
of Figure [2] for different values of the conditioned U3. Note that, due to the non 
continuity of Fu^Vz an d Fu 2 m 3 this copula is uniquely defined only outside the shaded 
areas. 
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5 Non continuous and non- simplified PCCs 




Figure 2 left panel: Scatterplot of the unspecified bivariate margin for a three- 
dimensional PCC of Cuadras-Auge copulas with parameters a 12 = 0.7 = 

"13|2, "23 = 0.3. 

right panel: Cdf of the bivariate conditional copulas of a trivariate MO 
copula with parameter A = 2 (see Appendix [U]). The cdf for U3 = 0.08 
is plotted solid (it is uniquely determined outside the light grey area), 
the dashed lines correspond to the cdf for w 3 = 0.28 (which is uniquely 
determined outside the dark grey area). 



Weakening the simplifying assumption 

Whereas the case of PCCs which are not absolutely continuous will rarely be relevant 
in practice, there is a wide range of applications for distributions involving non- 
simplified PCCs. 

Although the simplifying assumption sounds restrictive on the first glance, fairly 
general distributions can be obtained while keeping the simplified assumption for parts 
of the distribution. Staying in the parametric and absolutely continuous framework, 
the parameters of otherwise simplified PCCs can be made depending on the values 
of variables that are conditioned on. In cases where there is an underlying economic 
assumption for how covariates should influence the dependence, this can be done in 



the form of parametric models. Examples for this were studied by for example Patton 



(2006) and Bartram et al. (2007) , who considered models where the dependence 
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5 Non continuous and non- simplified PCCs 



parameter of the conditional copula depends on previous realizations of the time 
series, or Stober and Czado (2011) and Almeida and Czado (2011), who considered 



dependence of the parameter on an underlying Markov or AR(1) process, respectively. 



In particular, the Markov switching R-vine copula model of Stober and Czado (2011 
implies that the copula of a multivariate time series at each point in time is given by 
a discrete mixture of simplified PCCs, which will usually be of non-simplified form 
as for the mixture of Normals. 

When there is no a priori knowledge for how the dependence parameter should be 



influenced, non-parametric models as in Acar et al. (2011) can be applied. While this 



is fairly straightforward when conditioning on just one variable, it raises the question 
how interactions should be included when conditioning on multiple variables. For the 
elliptical and Archimedean distributions which we studied in Sections 3 and 4, we 
observe the following: 

• For elliptical distributions, the Kendall's r of conditional distributions does not 
depend on the values which are conditioned on. The distribution however can 
depend on these variables. Assuming zero means and correlations, and condi- 
tioning on realizations, Xk = Xk, ■ ■ ■ ,Xd = Xd the conditional distribution will 
only depend on a = Yl^krf- ^ the generator function gk(-) is monotonely 
decreasing (as for the multivariate Normal or Student-t distribution, and gen- 



erally if k > 2, see (Joe 1997, Section 4.9)), this implies in particular that the 
conditional distribution only depends on (^(x'^E^x/^). For this reason, one 
might consider to make the dependence parameter of conditional copulas de- 
pend on the likelihood of observations that are conditioned on, for data where 
the distribution appears to be close to the elliptical family. However, since the 
values of Kendall's r must not be affected — which are closely related to pa- 
rameter values for most well known bivariate parametric families — keeping the 
simplifying assumption will always be a close approximation in these cases. 

• For an Archimedean copula with generator function ip, we have observed that 
when conditioning on realizations Xk = Xk, ■ ■ ■ ,X^ = Xd the conditional dis- 
tribution will only depend on a = Yli=k < -P~ l { x i)- ^ n particular, this implies 
that it only depends on tp(a), which is the cdf of Xk, ■ ■ ■ , Xd, evaluated at 
Afc = Xk, ■ ■ ■ ,Xd = Xd- Thus, for data showing dependence behavior close to 
that of the Archimedean class, we recommend to consider analyzing dependence 
of the parameters of conditional copulas on this joint probability. 

As a simple example for how the conditional copula can depend on the values of 
conditioning variables in the Archimedean case, let us consider the three-dimensional 
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Frank copula. 

Example 5.1 (Three-dimensional Frank copula) 

Let U1..3 have cdf C Fra nk(uiS, a) = (p F (^(m; a) + ^ F l {u 2 ; a) + <fp X {uz\ a); a) with 
LT of the logarithmic series distribution: 



-a: 



Hog 



D + i 



a 



for a G (0, oo), 



i.e., the three-dimensional Frank copula. Then, the copula corresponding to the condi- 
tional distribution ofUi, U 2 | C/3 = u 3 is of the Ali- Mikhail- Haq family (LT of geometric 
distribution) 

C A mh(ui,u 2 ; (6(a,u 3 ))) = - — — — r 

1 -6>(a,w 3 )(l -u 2 ) 

with parameter 6(01,113) = 1 — e~ aU3 , thus depending on the conditioning value U3. 
How the strength of dependence between 1 , C/2 1 C/3 = U3 depends on the parameter a 
and the conditioning value U3 is illustrated in Figure^ 

Remark 5.2 (Extensions of Archimedean copulas) 

As for the multivariate Student-t copula, this decomposition allows for a natural exten- 
sion of the Frank copula by assuming different a parameters for the different bivariate 
copulas in the models, and similar extension arise for other Archimedean copulas. The 
possible range of dependence induced by this extension is illustrated in Figure [3| 





Figure 3 Density of the unspecified 13 margin of a trivariate meta-normal distri- 
bution with extended Frank copula with a%3 = 1 (r = 0.11), a 23 = 3 
(r = 0.21) and ai2|3 = 30 (right panel). For comparison, the left panel 
shows the density of a bivariate meta normal distribution where the de- 
pendence is characterized by a Frank copula with a = 3 (r = 0.21). 
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Figure 4 Kendall's r for the AMH-copula with parameter 9(a,u 3 ) = 1 — e~ aUA for 
different values of a and u 3 . 



6 Conclusion and summary 

While many popular statistical models are built by (sequential) conditioning, the 
implications on the properties of a distribution arising from this construction remained 
largely unknown. In addition, whether and under which assumptions popular classes 
of dependence models can be decomposed in this way was not investigated. 

In this paper, we filled in this gap by characterizing the most common classes of 
copulas in terms of their decomposability as a PCC. 

We showed that only the rf-dimensional Archimedean copula based on the gamma LT 
or its generalization into £ d can be decomposed using a PCC in which the building 
blocks are independent of the values that are conditioned on. For elliptical copu- 
las, the situation is more challenging. While the multivariate Normal and Student-t 
distribution are simplified PCCs, the conjecture that this is true for all elliptical dis- 
tributions does not hold. The Student-t distribution even is the only non-bounded 
elliptical distribution in which conditioning affects only the location and scale of the 
resulting conditional distribution but not the correlation matrix. Just as for the 
MTCJ in the Archimedean family, this distribution arises from a Gamma mixture for 
the square of the inverse scale parameter of a multivariate normal distribution. We 
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A Proofs involving Archimedean copulas - Theorem 



3.1 



have shown, that this makes it the only scale mixture of Normals which is a simplified 
PCC in dimension d > 4. 

Whereas the simplifying assumption for PCCs is convenient, it is often too restrictive, 
and also the assumption of dealing with absolutely continuous PCC is sometimes too 
strong. We illustrated this with several examples and demonstrated that in most 
applications, however, simplified PCC can be extended to adapt to the situation. 

While our classification results provide important insights into the properties of dis- 
tributions arising from conditioning, we believe that further research in this direction 
remains necessary, in particular considering how well general distributions can be 
approximated by simplified PCCs. 

Figure [T] does not include the extreme value copulas which is another general class. 
Extreme value copulas are used for data that are multivariate maxima or minima, and 
in this context, PCCs would not be considered. It can be shown that extreme value 
copulas which do not have independent subcomponents are not simplified PCCs; 
for example, given a trivariate extreme value copula C(iii, «2, which satisfies 
«2, w 3 ) = C t {ui 1 U2^u^) for all powers with t > 0, it can be shown that the cop- 
ula C*(v 1, V2', M3) obtained by conditioning on U3 = % has lower tail form v K, y(u3,v) 
as v — > 0, where K £ (1, 2) does not depend on M3 but the multiplying factor 7 is 
slowly varying in v and depends on u 3 . 



A Proofs involving Archimedean copulas - Theorem 3.1 



Because the proof is more intuitive, we first outline the case where (f is the LT of 
a positive random variable. In this case, there is a representation of the copula as 
a mixture of powers. The mixture distribution from which the Archimedean copula 



arises can be written as (c.f. Joe (1997, p. 86)): 

d 



H[G( Xj )r dF A [ 



a 



II [G(xj)] a f A (a)da, 



3=1 " u j=i 

where F A is the cdf of a positive rv, with corresponding density f A , and 



F(x) 



[G{x)] a dF A {a] = ^(-lnG(x)) 



is the common univariate cdf with ip A being the LT of A. Without loss of generality 
we can assume that F(x) = x on [0,1], then -Fi^xi^), Xj £ [0,1], j = l,...,d is 
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A Proofs involving Archimedean copulas - Theorem 



3.1 



a copula. Also G(x) = exp{—(p A (x)} on [0,1] is differentiable. With g = G', the 
marginal cdf of the last k variables and its density are: 

d 

\\ G{xi) f A (a)da, 



F(d~k+l):d(X(d-k+l):d) 
f(d-k+l):d(' x -(d-k+l):d) 



-k+1 



dx d - k +i ■■■dx c 

d 

n 



j=d-k+i 



9{ x j 
[G(x 3 



\d~k+l):d{^-{d-k+l):d) 
d 



a 



o 



i=d~k+i 



f A (a)da, 



and the conditional cdf of the first d — k variables given the last k is: 

d k F(x ll ...,x d ) 



^l:(<J-A)[(d-Jfc+l):d( x l: 



(d-k)\*(d-k+l):d) 



dx, 



d-kl 



■dx n 



/(d-fe+l):d( x (d-fc+l):d) 



/o°° n ti [G(x 3 T ■ « fc Utd- k+ i [G{xi)] a f A (a)d, 



a 



where 



Ia* (a) 



jo°° ak Ili=d-fc+i [G&T fA(a)da 

oo d—k 

Y[ [G(xj)} a ■ f A *(a)da, 

3=1 



k+1 [G( Xi )} a f A (a) ^ a k e aj:t =d _ k+1 HG {Xi) ) Ual 



f^nLw [OW f A (P)dfi 
In this case, A* has the same parameter form of density as A if 

f A (a; 77, 6) = e^V^/i^/C^, 9), (6) 

where h is a positive- valued function (it is not absorbed in the a 9 ~ l term only if it 
is a non-power function) and C(9,rj) is a (finite) normalizing constant. From above, 
A* has the same density form with parameters (r) — Ylt=d-k+i ^ n G(xi), 9 + k). The 
conditional copula does not depend on Xd-k+i, ■ ■ • > %d only if rj is a rate (or reciprocal 
scale) parameter, since Archimedean copulas are invariant to scale changes of the 
mixing distribution (Mai and Scherer 2012, p. 60). For rj to be a rate or inverse scale 
parameter of ^ only, h(a) must be a power of a. Hence f A is a gamma density, and 
F(xi : d) is a MTCJ copula. 

For the general case, where <p G Cd is not necessarily a LT, we prove the result by 
construction a functional equation. Let 

d 

c(ui, ...,u d ) = (ffy^^iuj 

3=1 
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A Proofs involving Archimedean copulas - Theorem 



3.1 



be an Archimedean copula. Let (Ui, . . . , Ud) be a random vector with this distribution. 
Suppose <p has support [0, So) where Sq = inf {s, (p(s) = 0} is infinite for a Laplace 
transform, but could be finite for (p G Cd which is not a Laplace transform. The case 
of finite support implies that (p(s) = for s > Sq. Let Ci...d_iw(ui, . . . , Ud-i\Ud) = 
dC(ui, . . . , Ud)/dud be the conditional distribution given Ud = Ud- We will show now 
that the copula for this is another Archimedean copula, say based on if), where if) G 
Cd-i- By differentiation, with h = —ip' and a = ip~ 1 (ud) G [0, Sq) with < Ud < 1, 

d 

Ci... d -i\d{ui, . . . ,u d -i\u d ) = /i^J^v 5 " 1 ^)) / H a )i 

0=1 

with jth (1 < j < d — 1) margin Fj(iij\ud) = h{(p~ l (uj) + a)/h(a) =: Note that 
h is monotone increasing by definition of Cd- Hence <y9 -1 (-Uj) = /i _1 (uj/i(a)) — a for 
1 < j < <i — 1 and the copula of the conditional distribution of Ui, . . . , Ud-i given 
U d = u d is: 

C*(v i, . . . , Vd-i, a) = h(^^ h^{v jh(a)) — (d — 2)aj j h{a) 

3=1 

Defining s = ip~~ l (v\ a) = h^ x {v h(a)) — a and v = ip(s; a) = h(s + a)/h(a), this is a 
Archimedean copula 

a) H h a); a) 

with generator function ^(-, a). As — )■ 1, a = (p~ 1 (ud) — > 0, and /i(0) = — y/(0) can 
be positive or infinite but not 0, by the definition of (p. 

Consider first the case where h(0) is finite and positive, and h'(0) is finite. The copula 
of the conditional distribution does not depend on Ud or a if and only if there is a 
continuous differentiable scale function 7(a) > such that 7(0) = 1 and 

if){s; a) = h(s + a)/h(a) = if>(s r y(a); 0) = h(s~f(a))/h(0); < s < s . 

Writing the above functional equation in h as h(s + a) h(0) = h(a) /i(s7(a)) and 
differentiating with respect to a yields 

h'(s + a) h(0) = h'(a) h(sj(a)) + h(a) h'(sj(a)) sj'(a). 

With a = if follows that 

h'(s) h{0) = h'{0) h{s) + h{0) h'{s) S7'(0), 
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A Proofs involving Archimedean copulas - Theorem 



3.1 



or 

ti(0)h(s) = h(0)[l-sj'(0)]ti(s) 

This has solution h(s) = h(0)[l - sY(0)] a where a = -ti(0)/[h(0)j'(0)]. Since 
h = —ip' must be decreasing, there are 2 possibilities 

(i) s = oo, 7'(0) < 0, a < 0, or 

(ii) 7 '(0) > 0, s = 1/Y(0), a > 0. 

By integrating h over s we obtain (p(s) = (1 — S7'(0))^_ +a , since ^(0) = 0. In case (i), 
we must have a < — 1 and in case (ii), I + a > d — lin order for ip £ £ d (see also 



Joe (1997, pp. pages 157-158)). In case (i) the obtained generating function support 



on [0, oo) and corresponds to the "standard" MTCJ copula, whereas case (ii), with 
bounded (p yields the extended MTCJ copula. 

If h'(0) or h(0) is infinite, the above is modified as follows. Let < e < s . There is 
a continuous differentiable scale function 7(a) > such that 7(e) = 1 and 

ip(s; a) = h(a + s)/h(a) = tp(sj(a); e) = /i(s7(a) + e)/h(e). 

Cross-multiplying and differentiating the above with respect to a, and then setting a 
to e yields 

h'{e) h{s + e) = h{e)[l - s-f'(e)]h'{s + e), < s < s - e. 

This has solution h(s + e) = h(e)[l — S7 ; (e)] a where a = — /i'(e)/[/i(e)7'(e)] so that 
ip(s; e) = h(s + e)/h(e) = [1 — S7'(e)] a . The conclusion is the same as above, because 
after integrating h to get if, one would conclude that this leads to (f'(0) and tp"{0) 
being finite. □ 

As an illustration for how the copulas of conditional distributions are derived when 
the conditioning is on more than one variable, let us consider the case where d > 4 
and 

/ d 
\j=i 

is an Archimedean copula. Let (U\, . . . , Ud) be a random vector with this distribution 
and let cij = ip~ l (uj) for j = d — 1 and j = d. Then, the conditional distribution of 
(U u C/ d _ 2 ) given U d _ x = u d -i, U d = u d is 

d 2 C(u 1 , . . .,u d )/du d ^u d 




Fl:(d-2)\d-l,d(Ui . . . , U d -2\U d -l, U d ) 



<f"(a d _i + a d )/[<p'(a d - 1 )<p'(a d )] 

¥" (j2jZi V 9 " 1 ^') + Od-l + a d 

^"(a^ + a d ) 
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B Proofs involving elliptical copulas 



The copula corresponding to this distribution is again Archimedean, based on a gen- 
erator ip. By differentiation, with h = ip", a = a^_i + a^, F 1 .( d _ 2 )\d-i,d{'\ u d-i-> u d) has 
jth (a < j < d — 2) margin Fj(uj\u d -i, Ud) = fo(</? _1 (/Uj) + a)/h(a) =: Vj. Note that 
h is monotone decreasing because ip is the generator of a rf-dimensional Archimedean 
copula. Hence, tp^iuj) = h~ l {vjh(a)) — a for 1 < j < d — 2 and the copula of the 
conditional distribution can be expressed as 

Ci:(d-2)|d-i,d(ui, . . . , Wd- 2 ; a) = /i h^ivjhia)) - (d - 3)aj / h(a). 

This is an Archimedean copula 4'(i J ~ 1 (vi 1 a)-\ \-ip~~ 1 (vd- 2 ', a)', a) with s = ip^iv; a) = 

h~ l {vh(a)) — a and v = ip(s;a) = h(s + a)/h(a). The same pattern extends to 
conditional distributions of Archimedean copulas with three or more conditioning 
variables. 



B Proofs involving elliptical copulas 
B.l Theorem 4.1 

In the remainder, we will use the following notation. A <i-dimensional Student-t 
distribution with mean vector 0, correlation matrix R and degrees of freedom v is 
denoted as td(0, R, v). Its pdf is f t ,d(-; R, v) and we write F t:d (-; R, v) for the cdf. 

Let us consider a d-dimensional random vector X = (X^,Xb) = (Xi,X 2 ,Xb), with 
A = {1,2}, distributed according to a multivariate Student-t distribution with v 
degrees of freedom, mean fi = (u-i, . . . , /i n ) T and scale matrix 



R = (Ri,j)i,j=i,...,d = ( ot R uj > where r ab = ( p 1B ) , Ra 



R\\ R\2 
R21 R22 



Let us define 
V41R : = 



R A \ B := diag(V A \ B ) 2 V A \ B diag(V A \ B ) 2 ,7<Xb) : = 



•1 + {1/v)k t b R- b 1 k b 



(v + d-2)/v 

then we have for the conditional distribution of X^ given X# = x#: 

P f , x „ ( -/xi|b(x b ) x 2 - // 2 |b(x b ) , \ 

V V^i|B -7(xb) v^V • 7(xb) / 
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B.2 Theorem 4.3 



c.f. Nikoloulopoulos et al. (2009, Lemma 2.2) Taking x 2 — > oo (x\ — > oo) yields 



j-i / | \ „ ( x i ~ ^i|b(x b ) 
Fi| B (xi|x B ) = F t ,i -j== ; rj f + a - 2 

\ v y i|s ■ T( x b) 
^ibO*^ xb) = F t ,i _ _ — - — -; v + d - 2 

so that we can now determine the corresponding copula: 

Cl,2|3:dOl,W2) = Cl,2|3:dOl,W2|xB) = i<12|3:d (-^s^fal I X b) > l X s) |Xfl 

= ^,2 (^.r'K; i/ + d - 2), F^r 1 ^; i> + d - 2); i? A |B, i> + d - 2) , 

where the additive constants and scaling factors cancel. This is a bivariate Student-t 
copula with v + d — 2 degrees of freedom and correlation matrix Ra\b and does not 
depend on x# anymore. □ 

B.2 Theorem 4.3 

The proof is similar to that for Archimedean copulas in that the same functional 
equation can be obtained. Without loss of generality, let us consider the case of a 
d-dimensional elliptical distribution with zero means and zero correlations. In this 
case, the density of the distribution is given by 

fl:d{ x i:d) = 9d{ x \ + ■ ■ ■ + x d ), with marginal density 
f(d-k+iy.d(*-(d-k+iy.d) = 9k(x 2 d -k+i + • • • + x d)- 

For fi.( d -k)\(d-k+iy.d{ x Hd-k)\ x (d-ky.d) we obtain 



/l:(d-fc)|(d-fc+l):d( x l:(d-fc)| x (d-/c):d) 



9d{x{ + ... + x d 

9k( x d-k+l + ■ ■ ■ + X d 



2^ 



For this distribution to be in the same location-scale family irrespective of the values 

of Xd-k+i, ■ ■ ■ ,Xd, we must have that for given J2i=d-k+i x i ^ Y,i=d-k+i x f tnere 
exists 7(x (d _ fc+1 ). d , x^_ fe+1) .J such that 

fl:{d-k)\(d-k+l):d{l ' X l:(d-fc) l X (d-fc+l):d) 00 fl:{d-k)\{d-k+l):d{^l:{d~k)\^{d-k+l):d)- 

With ^(d-k+iy.d = 0, a = Y? i= d-k+i x h t = EtN< and 6 ( a ) = 7 (x (tl _ fc 1 +1);d ,o) this 
implies that 

93 (t + a)=Z(a)-g 3 (j^Y (7) 
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B.3 Theorem 4.4 



where £(a) equals (71(a) times a constant depending on a and £(■), £(•) are differen- 
tiable scale functions. Since #3(0) = if and only if (73(a) = for all values of a, we 
must have (73(0) > 0. Using t — in Equation [7| we conclude that (73(0) is finite. 
Thus, we can define h(t) := |^ and obtain from Equation [7] that 

h(t + a) = a(a) ■ a 
Using that 5(0) = 1 by the definition of 5, 

- :V(t + „) = V(.).*(^) + *(.)-v(^).(^) 

a = : = a'(0) ■ fe(t) + h'(t) ■ (-t ■ 5'(0)). 

In other words, the function h must fulfill the differential equation 

(l + /3t)h'(t) = ah(t), 

where a = h'(0), /3 = 5'(0). From this, we obtain for (3 > that h(t) = (1 + /3t)f 
which corresponds to the elliptical generator of a Pearson type VII (scaled Student-t) 
distribution. For h to yield a well defined density in d dimensions, % must be given 
in the form s = — (V 4- d)/2, 1/ > 0, to ensure integrability with respect to t d / 2_1 . 

For < 0, a(t) = (1 — /3t)+, which leads to a well-defined density for § > —1 and is 
differentiable and thus a valid solution for % > 1. 

By integration, we obtain that the generator function for lower-dimensional margins 

9d-k{t) is proportional to (1 — /3t)+ 2 . Ergo, also the conditional distributions of the 
lower-dimensional margins remain within the same location-scale family for all values 
of the conditioning variables. □ 




B.3 Theorem 4.4 

A general scale mixture of Normals in dimension d can be written as (Xi, . . . , Xd) = 
(Zi, . . . , Zd)/VW where W is a random variable on (0, 00) with density fw, and 
(Zi, . . . , Zd) is multivariate Gaussian with zero mean vector and covariance matrix 
S. Without loss of generality, we can assume that all diagonal entries of S are 1, i.e., 
S is a correlation matrix. This implies for the d-variate generator gd that 

POO 

|E|- 1 / 2 y d (x'E- 1 x) = (27r)- rf / 2 |S|- 1 / 2 / w d/2 exv{-\wx!Y>- l x}f w (w)dw. 

Jo 
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B.3 Theorem 4.4 



Similarly, if H k is the leading k x k matrix of £, and X& = (x\, . . . ,Xk)', then by 
marginalizing out the last d — k components we get 



roc 

\£ k \- x l 2 g k {x! k Y.- l x k ) = (27r)- fe / 2 |S fc |- 1/2 / w k ' 2 exp{-^x' fc S fc l x k }f w {w) dw. 

Jo 

Hence, the univariate margin is 

POO 

f 1 (x 1 ) = 9l (xj) = (27T)- 1 / 2 / w 1/2 ex P {-lwxl}f w (w) dw. 



In the above, the density fw could be replaced with dFw in a Stieltjes integral. 
However, this would not affect the remainder of this proof; we omit it for notational 
convenience. Note that 

9l (0) = (2tt)- 1 / 2 E [W l '\ g d (0) = (2rc)- d / 2 E [W d ' 2 }. 

Let G be the cdf corresponding to the marginal generator g\. Then, the copula density 
for g d is 

c( Ul ,...,u d ;g d ,Z) = m- gi{{G -i {ui)) 2 ) ... gi{{G -i {ud)) 2 ) W 

with Xj = G~ l (uj). From this general form of the density, we will obtain two equations 
for the moments of the mixing variable W which will lead to necessary conditions on 
the conditional distributions of simplified PCCs in the elliptical class. Directly from 
(|sj> we get that 

This means that if W\, W2 are two different mixing variables, then the copula densities 
with fixed E are different unless the necessary condition of 

E [Wf /2 } E [W 2 d/2 ] 
E d [wl /2 } ~ E d [W* /2 } 
holds. For the second equation, let us consider 

c(0.5, . . . , 0.5, u- g d ) = \n^' 2 J^f) (10) 

where a > 1 is the (d,d) element of and x = Note that a is not a 

scaling factor of the marginal distribution but a function of the correlation matrix. 
In particular, we can obtain all a > 1 from correlation matrices 





d-2 








1 





1 - 
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B.3 Theorem 4.4 



where t d -2 is the (d — 2) x (d — 2) identity matrix. Using dx/du = l/gi(x 2 ), we 

2x 



determine the first derivative of (|10|) as 

isr 1 / 2 



ag' d (ax 2 



gaiax^g'^x 2 ^ 



Ygt^gxix 2 ) gf-\0)gKx 2 ) 

After taking the second derivative with respect to u and then setting x = (u = 0.5), 
all terms are except 



limlSl- 1 / 2 



ag' d {ax 2 ) _ ff d (ax 2 )ff / 1 (x 2 ) 

^(o)^^) ^woo 



g\(x 2 ) 



(2 7 r)- 1 |S|- 1 / 2 



a^(0) g d (0)gj(0) 
gt +2 (Q) 9^(0) 
aE (W^l 2 ) 



+ 



E (W d / 2 )E (I¥ 3 / 2 ) 
E d + 3 (WV2) 



Let us now consider the analogon for conditional densities. Let /i...d(x) be the density 
of (Xi, . . . ,Xd) and let /x be the density of X% or any Xj. Let 

/i...d-i|d(a?i, • • • ,Xd-i\x d ) = /i-d(x)//i(x d ). 

For this, we decompose x / S~ 1 x as (x*)'S 1 ~ 1 1 2 x* where x* = (x\ — (JidXd, ■ ■ ■ , %d-i — 
&d-x,dZd)' and En. 2 is the conditional covariance matrix of (Z\, . . . , Zd-i) given Z d - 
Writing the conditional densities in mixture form, 



'w) dw 



/i...d(x) = / w d/2 (p d (x.w;R) f w {w) dw 
Jo 

POO 

= (27r)- d / 2 |S|- 1 / 2 / /expf-IWE-^/p 
= (27r)- d / 2 |Sn. 2 |- 1 / 2 / w d / 2 exp{-i««*'Eii 2 x*}e^{-itw^}/ w (w)tZw J 



and 

= (27T)- 1 / 8 

so that 

fl-d-l\d( x l, ■ ■ > x d-lFd) — 



w 1 ^ 2 exp{ — ^wx d }f w (w) dw 



i(x d )(2n)-^/ 2 \^ 11 . 2 \- 1 ^ 2 / w^V 2 ex V {-\ W K*'Y;-l 2 K*}f w *(w)dw 

Jo 



'121 



is a scale mixture with mixing density fw*{w;x d ) = w 1 ^ 2 exp{ — ^wx 2 l }fw(w)/a(xd), 
where a(x d ) is a normalizing constant. We denote the random variable with this 
density by W*(x d ). 
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B.3 Theorem 4.4 



For d > 3, Equations ^ and (12) imply that a necessary condition for the copula 
corresponding to the distribution of Xi, . . . ,X d _i given X d = x d to be independent 
of x d is that 



J °° w d l 2 exp{-lwx 2 d }f w (w) dw- (j { 



00 w 1 / 2 exp{ — \wx 2 d }fw{w) dw 



d-2 



°°wexp{-iwi^}/iy(w) dw 



d-1 



(13) 



E[{W*(x d )} 



(d-l)/21 



is a constant over x d . Similarly, we obtain from (11) that 

aE(W*(x d ) (d+1)/2 ) 



\n- 1/2 



E^(w*{x d y/ 2 ) 



E{W*{x d )^l 2 )E{W*{x d fl 2 ) 

E d + 2 {w*{x d y/ 2 ) 



(14) 



must be equal to a constant /3(a). To rewrite these equations, let t = \x\ > and 
let V be a random variable with density fv(v) proportional to v 1 ^ 2 fw( v )- ^ Vt be 
random variable with density proportional to e~ tv fv{v); this is an Laplace transform 
tilt of the density of V with normalizing constant (pv{t), the LT of V at t. Note that 
Vt has finite positive integer moments for t > and Vq = V. 



Then, (13) can be rewritten as 



PCX 

Jo 



V 



(d-l)/2 , 



exp{— vt}fv(v) dv ■ ( exp{—vt}f v (v) dv 



d-2 



E[V t 



(d-l)/2. 



d-1 



( Jo°° t ' 1 ^ 2 ex P{ — v t}fv(v) dv 
which must be constant over t > while (|l4l) leads to 



E^[V t 1/2 } 



(15) 



aE(^ (d+1)/2 ) E(v/ d - 1)/2 )E(V; 3/2 ) 



E d +!(^ 



1/2, 



+ 



Ed+2(y t 1/2 ) 



(16) 



which, for all t > 0, must be equal to a constant /3(a). This implies the following 
recursive relation for the moments of V t : if we know that (15) is constant for d = k 



and d = 4, then it is also constant for d — k + 2. Thus, it is sufficient to show that 



(15) is constant for d = 3 and d = 4, or d = 3 and d = 5. 



Let us now consider case a), where the copula C is a simplified PCC in d > 4. From the 
three-dimensional marginal distributions, we obtain that (15) is constant for d = 3. 
By conditioning on X 4 = x^,X 3 = x 3 , we conclude with a similar calculation as for 



(12) that 



E[V 2 } ■ E 2 [V t 1/2 } 
E*[V t ] 



const. 
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C Trivariate Marshall-Olkin (MO) copula 



This, together with (15) being constant for d = 3 implies that (15) is constant for 
d — 5 and thus for all d > 3. 

In case b), where the copula is a simplified PCC for all E, we know that for a three- 
dimensional marginal distribution, (16) holds for all a > 0: 

aE i{W*(x 3 )y/ 2 } E[{V^(x 3 )} 2/2 ]E[{V7*(x3)} 3 / 2 ] , 



E 4 [W* ' (x 3 )] 



E 5 [W* 1/2 (x 3 )] 



Thus, E [V 2 ]/E 4 [V t 1/2 ] and E [Vt]E [V t 3/2 ]/E 5 [V t 1/2 } must be constants over t. Together 
with ((TBI) being constant in t for d = 3, this means that (15) is constant for d = 4 



and thus for all d > 3. Note that, more precisely, we only require two different values 
of a in (16) for the argument above. 

Summing up, we obtain that with m = (d — l)/2 the moments of V t are connected to 
the moments of V via 

E[V t m ] = a m (t)E[V m ] 

for alU> 0, m= 1,1.5,2,2,5,..., where a(t) = E 2 [V t 1/2 }/E 2 [V 1 / 2 }. 

With all of the positive integer moments of V and Vt existing, the Laplace transforms 
of V and Vt, for < s < St, where the constant St may depend on t, can be written 

as 

oo 

< Pv (8) = i + Y^(-iYE\y i ]8 i /ii 1 



1=1 



1 + j^(-l)'E [V;VA! = 1 + X](-1) 4 E [V^a^/il 



i=i 



i=l 



Hence v^Vtl 5 ) = V 9 v( sa (^)) i n a neighborhood of for the Taylor series expansion 
of the LTs about 0. By Feller (1971, Section VII. 6), the Taylor series in a positive 
neighborhood of uniquely determines the distribution. Hence Vt = a(t)V for t > 0. 
The combination that Laplace transform tilting of the density leads to a scale-changed 
random variable, implies that V has a gamma density (Marshall and Olkin 2007 p. 
576, Theorem 18. B. 6). Hence, also W is Gamma distributed, and the corresponding 
scale mixture is the multivariate t-distribution. □ 



C Trivariate Marshall-Olkin (MO) copula 

To determine the bivariate conditional distributions of a three-dimensional MO cop- 
ula, we work with the parameterization of Mai and Scherer (2012, Chapter 13)] The 
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C Trivariate Marshall-Olkin (MO) copula 



three-dimensional MO copula is the survival copula of the rvs defined as 

Xi : = min |-Ei, £12, E i3 , -E123 j 
X 2 := min j_E 2 , E i2 , E 23 , E 123 ^ 
X 3 := min j_E 3 , E 13 , E 23 , E 123 ^, 

where Ej, I C {1, 2, 3} are independent and exponentially distributed. For simplicity 
let us assume that all rate parameters are equal, i.e., A/ = A, V/ C {1,2,3}. This 
implies that P(X 3 = Ej\X 3 = x 3 ) = 1/4 for all Ei,I C {1,2,3}, and thus we can 
determine the conditional survival distribution of Xi,X 2 \X 3 = x 3 as 

F(x 1 ,x 2 \x 3 ) = P(X l > Xi,X 2 > x 2 \X 3 = x 3 ) 

= P ( X i > x *> X * > x ^ = x ^ X 3 = Ej) ■ P(X 3 = Ej\X 3 = x 3 ) 

3e/|/C{l,2,3} 

= \- P(X 1 >x 1 ,X 2 >x 2 \X 3 = x 3 ,X 3 = E I ). 

3e/|/c{l,2,3} 

If X 3 = x 3 and X 3 = Ei 23 , we know that x\, x 2 > E 123 > E 23 , E\ 3 . Using this, 

P(Xi > x u X 2 > x 2 \X 3 = x 3 ,X 3 = E 123 ) 

= t xl , X2 < X3 P(E 1 > x u E 2 > x 2 ,E 12 > max(xi,x 2 )), 

since X 2 > x 2 implies E 2 , E 12 > x 2 and X 1 > X\ implies E ± , E 12 > xi, which further 
yields E i2 > max(xi, x 2 ). 

Let X{ X >y} be the random variable with cdf P(X < x\X > y), then we obtain 

P(X 1 >x±,X 2 > x 2 \X 3 = x 3 ,X 3 = E 23 ) 

= t X2 < X3 P(E 1 > x 1 ,E 2 > x 2 ,E 12 > max(xi,x 2 ),E 123{El23>xs} > x ± , 

El3{E 13 >x s } > Xl), 

P{X 1 > Xl ,X 2 > x 2 \X 3 = x 3 ,X 3 = E 3 ) 

= P(E 1 > xi,E 2 > x 2 ,E 12 > max(x 1 ,x 2 ),E 2 3{E 2 3>x3} > x 2, 
Ei 3 {e 13>Xs } > xi, E 123 { El23>X3 y > max(xi,x 2 )), 
P(X 1 >x ± ,X 2 > x 2 \X 3 = x 3 ,X 3 = E 13 ) 

= l Xl < x . a P(E 1 > X!,E 2 > x 2 ,E 12 > max(xi,x 2 ),E 123{El23>xs} > x 2 , 

E23{E 23 >x 3 } > X 2 ). 
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Putting together the conditional probabilities, the conditional survival function is 

F m (x 1 ,x 2 \x 3 ) = I e -A(*i+*a+™«(*i.*2)) 

1 4 i 71 -A(xi+X2+max(xi,a;2)-3a;3\ 

^■Xl,X2<X3 ^ ' ^■X1,X2>X3 C ) 
"i" lrri<a;3,X2>a;3 • 2 • e ^ 3 ^ 4" lrr 2 <2:3,2:i>:r3 ' 2 • e ^ 3 ^ . 

For the generalized inverse of the conditional survival function of e.g. = x 3 , 

this implies that 

Given this, we can evaluate the copula of the bivariate conditional distribution of 
Xi,X 2 \X 3 = x 3 to see that it depends on the value of x 3 in the areas where it 
is uniquely determined. To arrive at the conditional copula of the trivariate MO 
copula, and not of the corresponding distribution, we further have to transform x 3 . 
For the marginal survival function we obtain that F 3 (x 3 ) = e~ 4Xxs , with inverse 
F 3 - 1 (v) = -\og(v)/(4\). 
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